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In this paper we analyse the possibility of constructing singularity-free inhomogeneous 
cosmological models with a pure radiation field as matter content. It is shown that the 
conditions for regularity are very easy to implement and therefore there is a huge number 
of such spacetimes. 
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1. Introduction 

Singularity-free cosmological models were quite popular until the sixties. At that 
time there were attempts to circumvent initial and final singularities, for instance, 
by introducing additional fields, such as Hoyle and Narlikar's C field. Such models 
had the drawback of requiring continuous creation of matter U4 The interest for 
that line of research decayed in the seventies when most of the famous singularity 
theorems were published (cfr. for instance 1^,1^. Under very general and physically 
realistic conditions these theorems seemed to push cosmological models towards the 
appearance of singularities, that could be of the Big Bang kind or not. Although 
most quantum fields, such as Yang-Mills, do not follow them, there was nothing to 
object to energy conditions. Less criticism would receive causality conditions, which 
require the absence of closed causal curves. Or the generic condition, which roughly 
states that the universe is not empty El. 

However, in addition to those general and physically acceptable requirements, 



February 7, 2008 2:47 WSPC/INSTRUCTION FILE 



radiationl 



2 



there were other conditions that might be overcome. These were the existence of 
some causahy trapped set in the spacetime, such as closed trapped surfaces El or 
an achronal set without edge The fulfillment of the latter conditions is not so 
obvious as for the former ones. 

But it was not until the nineties when the first singularity-free cosmological 
models began to appear El These models were not realistic, since they were in- 
homogeneous spacetimes with cylindrical symmetry, but they pointed out that on 
diminishing the symmetry non-singular behaviour might appear. Since these models 
were exact solutions in a framework where they are fairly difhcult to obtain, it is not 
strange that the list of such non-singular models was rather short until the end of 
the centurylZl. In fact these spacetimes belonged to larger families among which the 
absence of singularities required a certain amount of fine tuning El. Consequently it 
was tempting to consider singularity-free cosmological models as a curiosity. 

However a large family of regular cosmological models was recently found ^ 
and did not depend on a finite number of parameters but on two nearly arbitrary 
functions. This result could mean that although our universe has an initial singu- 
larity, it is yet to be explained why alternative, and abundant, possibilities were not 
considered. 

The aim of this paper dates back to Senovilla's original solution. This spacetime 
had the realistic property of being filled with an incoherent radiation perfect fluid, 
which is the standard mean matter content for a universe in its early epoch. It is 
intriguing that when the Senovilla spacetime was extended to include other models 
with also a linear equation of state again radiation fluids were the only ones 
which remained non-sing ularlini. What we plan to discuss now is the role of radiation 
in non-singular behaviour. But instead of considering incoherent radiation, we shall 
focus on pure radiation. 

With this aim in mind we devote the next section of this paper to the Ein- 
stein equations for radiation fields in inhomogenous cosmological models with a 
two-dimensional group of isometrics acting on spacelike orbits. Section 3 comprises 
three ways of simplification of the equations. One of them, corresponding to models 
with a transitivity surface element with lightlike gradient, incorporates additional 
symmetries to the spacetimes. In Section 4 the analysis of the possible singularities 
of the families described in the previous sections is completed. A discussion of the 
results is included at the end of the paper. 

2. Inhomogeneous radiation cosmological models 

We may describe pure radiation by its conserved energy-momentum tensor. 



where A: is a null vector field describing the lightlike 4-velocity of the radiation, 
k^kf^ — 0, and $ is a positive function related, but not equal, to the radiation energy, 



= ^kt'k", 



(1) 



I 
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since this function could be reabsorbed by the vector field k without disturbing the 
properties of the tensor. Bear in mind that radiation fields are characterised by 
having all energy-momentum tensor invariants equal to zero. 

The energy-momentum conservation equations may be written, 



T';^ = <i>„k''ki' + + $efc'" = 0, (2) 

in terms of the expansion Q — fc.^j, and acceleration, — fc^^fc" of the lightlike 
congruence. Since the acceleration is orthogonal to the radiation velocity, a might 
have a component parallel to k and a spatial component orthogonal to k. The 
existence of the latter component would mean that the equations would split into 
components parallel and non-parallel to k and $ would be equal to zero. Since this 
is not our case, a is proportional to k and |(2J) is a single-component equation. 

Now we turn to the geometric description of the spacetime. Considering that 
except for some spherically symmetric spacetimes ' every other singularity-free cos- 
mological model has been found within the Ansatz of an Abelian two-dimensional 
orthogonally transitive group of isometries acting on spacelike orbits, we shall im- 
pose that our models share these symmetries. We choose as generators for the 
symmetries two commuting fields, that we write as {dy,dz}, having in mind the 
coordinates that we want to use. We take them as mutually orthogonal. There is 
much freedom for choosing the rest of the coordinates, which we shall name i, x. 
Since every surface admits and isothermal parametrization for which the metric is 
conformally flat, we take as our t and x the coordinates in such parametrization. 
With this choice of coordinates, the metric can be written in a suitable form^^, 

ds^ = e^^'i-dt^ + dx") + p^e^^'dy^ + e-^^'dz^, (3) 

depending on three functions, K, U and p, of just t and x. 

As we shall see, there is no possibility for having an axis in these models. There- 
fore we reject the interpretation of the isometries as cylindrical symmetry and there 
is no reason a priori for restricting the ranges of the coordinates, 

— cxD < t,x,y, z < oo. (4) 

As an Ansatz we take a velocity of the radiation with no projection on the 
Killing fields, that is, 



fc± = e-^(at±9,), (5) 

where the plus (minus) sign stands for growing (decreasing) coordinate x. The 
normalization factor has been chosen so that fc_) = —2. 

For such velocities the energy-momentum conservation (j2l is written as 

$t±$, + $(^^^ + 2(i^t±if,)l =0, (6) 
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which allows a rather simple form in terms of ingoing or outgoing lightlike coordi- 
nates, 

(at±9,){$e2^^p} =0, (7) 

for which integration is straightforward, 

'S>e^^p^h(tTx), (8) 

in terms of a free function h of one outgoing (ingoing) lightlike coordinate. 

Finally, the set of Einstein equations for pure radiation can be cast in the form 

Utt-U^x + ^{UtPt~U^Px) = 0, (9a) 

Ptt - Pxx = 0, (9b) 
Ktpx + Kxpt = Ptx + Utpx + UxPt + 2pUtUx T h, (9c) 
Ktpt + KxPx = ^^i±^ + UtPt + axPx + P {u! + Ul) + K (9d) 

Kxx - Ku + ^^^^ ~ + Ul Vl = 0, (9e) 
P 

where the integrability condition for (|9cl9dl) is the Euler equation |(7|l and the 1- 
D inhomogeneous wave equation for iiT is a consequence of the rest of the 
differential system and therefore can be dropped. The remaining set of differential 
equations for pure radiation models is then reduced to a 2-D wave equation for [/, 
a 1-D wave equation for p and a quadrature defining K . 

3. Simplification of the equations 

The system of differential equations may be simplified by taking p as a coordinate 
for the equations. Three possibilities arise according to the character of the gradient 
of p: 

• If the gradient is lightlike, pt = spx, = 1, the 1-D wave equation is 
trivially satisfied, pxx = £ptx = spxt = Ptt- We have then p = p(t + ex) and 
the equations defining the quadrature, 

(Kt + ei^x)- = ^ + {Ut+ eUx)- + 2eUtUx T e^e^^, (10a) 
P P P 

{Kt+eK,)^^^ + {Ut+eUx)^ + U? + C/J + c&e^^' , (10b) 
P P P 

allow us to draw rather constraining consequences, 

$e2^(l±e) + (C/t-eC/,)2 = 0, (11a) 

KtTK^ = — + m + '2p^+p , (lib) 

Pt Pt Pt 

since e = =f1 in order to have $ 7^ and the gradient of U is also lightlike, 

Ut — tUx, and parallel to the gradient of p, but not to the velocity of the 

radiation. 
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No further restrictions arc implied for U, since it satisfies Utt = Uxx and 
the remaining condition, 

= PxUx - PtUt^ (gradJJ, gradp) (12) 

is trivial. 

For K the 1-D wave equation H9e|) is no longer inhomogeneous, 

Ktt - Kxx = 0, (13) 

and we may write the solution in terms of ingoing and outgoing lightlike 
coordinates, u — x, v — t ± x, 

K{u,v)^f{u)+g{v), (14) 

but just / is constrained by the quadrature, 

where the prima denotes derivation with respect to u. 

In principle, we have four independent functions, /, g, p and U and the 

system is solved by prescribing 



h = 2p'{f' -U')- p" -2pU'^. (16) 
The metric can be therefore written as 



ds2 ^ -e^f+^9dudv + p^e^^dy^ + e'^^dz^, (17) 
in terms of lightlike coordinates. 

But we see that / and g appear due to just coordinate choices. We may 
eliminate them by respectively redefining u and w as J du e^^ and J dv e^^ 
and writing the metric. 



ds^ = -2dudv + F^dy^ + G^dz^, (18) 

in terms of just two independent functions, F and G, of u. 
In this simpler frame the function $, 

just imposes the positiveness requirement. 

It is clear that these spacetimes possess more isometrics than those which 
were required from the beginning. Besides the obvious Killing fields, 
{dv,dy,dz} there are two additional independent isometry generators, 
{ydv + J F~'^ dy, zdv + J G^^ dz} and therefore the isometry group is a 
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6*5 acting on null orbits. The spacetimes are not however homogeneous, 
since the action of the group of isonictries is not transitive. The Petrov 
type is generically I. 

The Killing field d„ is covariantly constant and therefore these models are 
plane-fronted gravitational waves with parallel rays spacetimes HSI, El, El. 
There are examples where the isometry group is even larger, for instance 
when F = G, there is an additional independent Killing field, ydz — zdy. In 
this case the spacetime is conformally flat. 
The results may be extended to non-diagonal metrics. 



ds^ = -2dudv + F^dy^ + 2Hdydz + G^dz^, (20) 

by just introducing another function, H , of u. 
• If the gradient of p is spacelike, we may choose p as a coordinate, p = x. 
Equation (j9b|) allows us to keep the isothermal parametrization. The system 
thereby becomes much simpler. 



Utt - C/.. - — = 0, (21a) 

X 

Kt = Ut + 2xUtU^ T h, (21b) 
K^^U, + x (C/2 + Ul) + h, (21c) 

which tell us that radiation fields can be generated from vaccuum space- 
times by adding a term ^H{t =F x) to the conformal factor K . 
This is the only case where a true axis might appear at x = 0. But on 
applying the flatness condition in the vicinity of the axis El 



lim (g^adA,gradA) ^^,(,_^)^^^^^^^ A = (a., 9.) = xV^, (22) 

a;— >0 4Z\ 

assuming that the metric functions are regular at a; = 0, the quadrature 
equations for I|21bl21c|l . would imply 



Kt{t,Q) = Ut{mTHt[t), K,{t,Q) = U^[t,0)TH,{t) (23) 

naming h = H'. Therefore K(t,0) = U{t,0) + H{t) and condition ^ 
cannot be satisfied for all times unless H is constant. 

Furthermore, this means a problem in our parametrization, since the func- 
tion h{t =p x) must be positive (negative) for positive (negative) values of 
x. Therefore, if h(tQ =F xq) is positive for some values xq, to, it should be 
positive for all values of t and h would be everywhere positive. It is obvious 
that there is no place at the same time for positive and negative values of x 
and we have to resort, although there is no axis in strict sense, to the axial 
interpretation of the models and admit that x ranges from zero to oo and 
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it is a sort of radius. The coordinate y would be interpreted as as an angle. 
In such case, Ht = H' = h must be positive everywhere and therefore the 
conformal factor H is a growing function of its only variable. 
• Similarly, if the gradient of p is timelike, we may choose p as a coordinate, 
p — t without losing generality, 

Utt - + y - 0, (24a) 
K, = U, + 2tUtU, T h, (24b) 
Kt^Ut + t {Uf + Ul) + h, (24c) 

for which the results of the previous case are also valid. 



4. Geodesic completeness of the models 

The most usual definition of absence of singularities is causal geodesic completeness. 
This means that every timelike and lightlike geodesic, that is every physical trajec- 
tory of a test particle in the spacetime in the absence of external non-gravitatory 
forces, can be extended from — oo to oo in their affine parametrization, = x'^{t), 



+ T'^pi'^xP = 0, g^^x'^xP = -(5, -oo < t < cx) (25) 

where 5 takes the value zero for lightlike and one for timelike geodesies and the 
dot stands for derivation with respect to the affine parameter t. It is called affine 
parameter because it is unique but for affine transformations. 

Causal geodesic completeness implies a simpler definition, regularity of the cur- 
vature scalars, but it is not true the other way round. Taub-NUT spacetime ^1 
is a nice counterexample for this Moreover, spacetimcs where curvature scalars 
are regular may be incomplete because they can be extended to larger manifolds. 
Therefore in order to check geodesic completeness we are to analyse geodesic equa- 
tions. 

We have already seen that our models are classified in three sets according to 
the character of the gradient of p: 

• Lightlike gradient of p: In this case the analysis of the geodesies becomes 
quite simple, since there are extra isometrics, which allow us to reduce 
the order of the geodesic equations. We just need the three commuting 
generators {dy, dy,dz} to derive independent constants of geodesic motion, 



V = 


{dy , X) 


— Qvfi^^ 


= -u, 


(26) 


Y = 


{dy,x) 






(27) 


Z = 






= G^i, 


(28) 
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The first equation simply implies that u is an acceptable affine parameter 
for the geodesies, 

u = uo-Vt. (29) 
The second equation can be integrated for all values of t or u, 

y^yo-^J duF~\u). (30) 

provided that is integrable for all values of u. 

The third equation is also integrable if is integrable for all values of u, 

z ^ zo- ^ J duG^^iu). (31) 

The remaining equation of geodesic motion needs not be derived from the 
Christoffel symbols but from the condition of affine parametrization, 

S = -2Vv -Y'^F-'^ - Z^G-^, (32) 
from which we do not get any additional restriction, 

^ = ^{S + y'F-' + Z'G-^}- (33) 

since v is integrable provided that y and z are integrable. 
These results may be summarized in the following theorem; 

Theorem 1. An spacetime with a metric of the form H18|) filled with a pure 
radiation field is causally geodesically complete if and only if the metric 
functions g^^ and g^^ are integrable for all values of the light coordinate u. 

Examples of such spacetimes are easy to provide. For instance, if we resort 
to powers of the coordinate u, F{u) — \u\°'^^, G{u) — \u\^^^, the regularity 
requirement is fulfilled ii a, /3 < 1. Since $ has to be positive, an additional 
condition is imposed, < a, /3 < 2. This example is intriguing since it states 
that half of the models are regular whereas the other half are singular. 
Although the previous reasoning has been made for specific functions, it 
describes the behaviour of general functions F and G on becoming zero. 
In fact, since equation 1)19(1 imposes that the graphic of at least one of the 
functions must be convex, it is clear that the graphic of this function is to 
cut the axis at one point at least. 

This poses a serious drawback on these spacetimes, since the previous anal- 
ysis points out that the metric functions F^ and G^ are continuous, but 
not differentiable. We are dealing with topological manifolds instead of dif- 
ferentiable manifolds. However, metric invariants do not see this fact, since 
they are all null. 
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• The analysis of models with non-lightlike gradient of p is quite more compli- 
cated, since generically we do not have additional symmetries. Fortunately, 

we may resort to some theorems which have been explicitly devised for G2 
1 o 

spacetimes , 

Theorem: A diagonal Abelian orthogonally transitive spacetime with 
spacelike orbits endowed with a metric in the form (PJ with metric 
functions K, U, p, where p has a spacelike gradient, is future causally 
geodesically complete provided that along causal geodesies: 

(1) For large values of t and increasing x, 

(a) {K~U -\np)t + {K~U - \np)^ > 0, and either {K - U - 
Inp), > or \{K-U -In p),\ < {K~U-\iip)t + {K-U-liip),. 

(b) Kt + Kx> 0, and either K^>0 ot \K^\ < Kt + K^. 

(c) {K + U)t + {K + U)^ > 0, and either (K + U)^>Oot\{K + 
UU < {K + U)t + {K + U)^. 

(2) For large values of t, a constant b exists such that 

K{t,x)-U{t,x) 1 

2K{t,x) I > - In |t| -t- 5. 

K{t,x) + U{t,x) + In p{t,x) J 

For geodesies pointing to the past the theorem imposes the same conditions 
but for a sign in the time derivatives. 

We just have to check the restrictions on the metric functions that such 
conditions require: 

(1) Since h is positive for radiation fields and positive values of x, we may 
reach some conclusions: 

(a) [K -U - \np)t + {K -U -\n p)^ > x {Ut + U^f - 1/a; is positive 
unless the derivatives of U decrease too quickly. Such possibility is 
avoided if we impose 



x^-'\U^ + Ut\^0, (34) 
for large values of t and x. 

On the other hand {K — U — Inp)^ > x {U"^ + C/J) — \/x \s also 
positive under the same conditions. 

(b) Kt+Kx > Ut+Ux+x {Ut + Ux) is again positive for large values oft 
and X if is satisfied, as it happens with > Ux + x {U^ + U"^) ■ 

(c) The same restrictions are valid for {K + U)t + {K + U)x and \{K + 
U)x\- Therefore l)34|l suffices for the satisfaction of the first part of 
the theorem. 

For past-pointing geodesies, as it has already been told, condition 1 is 
the same for large negative values of t, but for the sign of the time 
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derivatives. Therefore, the corresponding requirement for past-pointing 
geodesies is 

x^-'\U^-Ut\^Q, (35) 

for small values of t and large values of x. 
(2) On a; = we aheady know that K{t, 0) - C/(i, 0) = TH{t) and H is a. 
growing function. Therefore we may encounter a problem for outgoing 
radiation unless H grows slower than a logarithm for large times, 

TH{t)>-\n\t\+b. (36) 

Out of a; = the first condition of this set is satisfied, 

K{t,x)-U{t,x)^TH{tTx)+ f x' {Utit,x)^ + U^it,x)^} . (37) 

Jo 

since we are adding just a positive term to the previous case. 

For the second condition on a; = 0, we obtain another constraint for 

large values of t, 

K{t, 0) = U{t, 0) T H{t) > In \t\ + b, (38) 

and out of .t = it will be also satisfied. 
Finally, for the third condition we may require, 

2U{t,0)T H{t) > -lii\t\+b, (39) 

always for large values of t. 

These conditions on outgoing radiation may be summarized in the following 
way for large values of i. 



H{t) < ln|t| +6 
U{t,0)>b-^\n\t\+H{t), (40) 
since H grows obviously faster than H/2. 

For ingoing radiation condition H36|) is always fulfilled, since if is a growing 
function. The only requirement which remains is 

Uit,0)>b~lln\t\-^, (41) 

for large values of t. 

Similar restrictions are imposed for past-pointing geodesies for large nega- 
tive values of t. The roles of outgoing and ingoing radiation are obviously 
exchanged. 
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For outgoing radiation the only restriction in order to achieve causal 
geodesic completeness to the past is 

Uit,0)>b-Un\t\ + ^^, (42) 

for small values of t, since H/2 decreases to the past more slowly than H. 
On the other hand, the corresponding condition for ingoing radiation is 

H{t) > -ln|t| +& 
U{t,0)>b-^\n\t\-H{t), (43) 

for small values of t. 

Summarizing the previous results, we may enclose them in two statements: 

Theorem 2. A diagonal Abelian orthogonally transitive spacetime with 
spacelike orbits with a metric in the form Q and a spatial gradient of p 
and a pure outgoing radiation field as matter content is causally geodesically 
complete if: 

(1) x^~'^\Ux ± [/f I 7^ for large values of \t\ and x. 

(2) H{t) < In \t\ + b and U{t, 0) > - i In |t| + H(t) for large values of t. 

(3) U{t, 0) > 6 - i In |i| + H{t)/2 for small values of t. 

Theorem 3. A diagonal Abelian orthogonally transitive spacetime with 
spacelike orbits with a metric in the form iPJ and a spatial gradient of p 
and a pure ingoing radiation field as matter content is causally geodesically 
complete if: 

(1) x^''^\Ux ± C/f I 7^ for large values of \t\ and x. 

(2) U{t, 0) > 6 - i In |t| - H{t)/2 for large values of t. 

(3) H{t) >-ln\t\+b and U{t, 0) > & - i In |t| - H{t) for small values of t. 

The meaning of such results seems clear. There must be an upper limit for 
H in the future for outgoing radiation in order to prevent a huge amount of 
energy that could form a singularity. A symmetric reasoning can be applied 
to ingoing radiation. 

H can be any growing function, but U is restricted by the 2-D wave equation 
(|21a|) . Following we may write the general solution of this equation in 
terms of initial data, U{0,x) ~ f{x), Ut{0,x) = g{x), 

U{t,x) 



where v 



1 

2^ 



dT 



tg{v) + f{v)+tf{v[ 



It' 



XT COS ( 



(44) 



= Jx^ + <2 



2xtT COS ( 
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Since we are interested in the values of C/(i, 0), the expression we need is 
much simpler, 



which allows easy implementation of the regularity conditions. 
For instance, we can construct non-singular radiation models just prescrib- 
ing a function U which grows for large values of |<| and x and a function H 
that does not grow or decrease too fast, that is, bounded by a logarithm. 
A simple example of a way to comply with the regularity requirements 
for U is shown in 1^ for polynomial initial data: if f{x) — '^aiX^ and 
g{x) = ^ bjX^ are polynomials in x of degrees m, n, respectively, U fulfills 
the regularity conditions if am > and m > n -f 1 or if m = n + 1 and 
{m + l/2)a,n >\bn\. 
• We shall not consider causal geodesic completeness in the case of a time- 
like gradient of p since the logarithm term of the time coordinate in the 
regularity theorems would mean a singularity at t = 0. 

5. Discussion 

We have discussed in this paper the possibility of constructing regular inhomogenous 
spacetimes with a pure radiation field as matter content. It has been shown that 
such cosmological models are very easy to construct in terms of a nearly arbitrary 
growing function H, which is bounded by a logarithmic function, and a solution of 
a wave equation, U, with rather few limitations. The conclusions are close to the 
results for stiff perfect fluids: pure radiation inhomogeneous cosmological models 
are an inflnite family and are not negligible in the set of much models. 

Pure radiation fields satisfy all energy conditions and these models are causally 
stable since the function t has a positive gradient in all the spacetimeEl Therefore 
these models fulfill weaker causality conditions suchs as the chronology conditions. 
As it happened in most regular cosmological models, the way to elude singularity 
theorems is the lack of trapped sets. 

The existence of such an amount of non-singular cosmological models fulfilling 
all energy and chronology conditions suggests that they are more abundant than it 
was thought before. If they had not been found previously it was due merely to the 
techniques used to generate them, not to their scarcity among inhomogeneous cos- 
mological models. It remains however an open question whether such non-singular 
behaviour is just a peculiarity of G2 inhomogeneous models and if these results can 
be extended to smaller groups of isometries or to other matter contents. 

As a curiosity, it is interesting to notice that perfectly regular spacetimes can 
be turned singular by introducing radiation. For instance, the metric, 




(45) 



ds^ = e'^-* {-dt^ +dr^}+ r'^dcf^ + dz^ , 



(46) 
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can be generated from Minkowski spacetime by adding a conformal factor. It is easy 
to check that axial geodesies are not complete, 

r = 0, 4) = 0, S = e-H'^-Z'^, (47) 



i = VS + Z2e* ^ t(r) = - In (V^ + Z2(fc - r)) . (48) 
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